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An  analytical  treatment  of  the  “strong*  interaction  problem 
in  hypersonic  viscous  flow  on  a  noninsulated  flat  plate  is  pre¬ 
sented,  using  the  method  of  “similar*  solutions  of  the  compres¬ 
sible  boundary  layer  equations.  Recent  experimental  data 
which  confirm  some  of  the  theoretical  results  are  also  discussed. 
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«  coefficient  of  viscosity 

coefficient  of  heat  conductivity 
specific  heat  at  constant  pressure 
gas  constant 
total  energy 
ratio  of  specific  heats 
defined  in  Eq  5 
defined  in  Eq  5 

the  constant  parameter  in  Eq  1 
similarity  variable 
boundary  layer  thickness 
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A,  B,  C,  D 

tJ 


local  skin  friction  coefficient 


local  heat  transfer  coefficient 

free  stream  Reynolds  number  on  the 
basis  of  x 

defined  in  Eqs  16-19 
exponent  m  the  viscosity  law 


Subscripts  “o',  “1",  *«,*,  *w”  refer, 
respectively,  to  the  stagnation  condi¬ 
tions,  the  conditions  at  the  edge  of  the 
boundary  layer,  the  conditions  in  the 
free  stream,  and  the  conditions  on  the 
surface. 
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Introduction 

The  study  of  hypersonic  viscous  gas  flow  past  a  flat  plate 
constitutes  one  of  the  best  sources  of  information  on  the  nature 
of  hypersonic  interaction  between  the  viscous  and  inviscid  ef¬ 
fects.  Serious  efforts,  both  analytical  and  experimental,  have 
recently  been  directed  towards  the  understanding  of  this  in¬ 
teresting  problem.  Shen  [1]  (Numbers  in  brackets  refer  to  the 
bibliography  at  the  end  of  the  paper.)  first  attempted  in  1949 
an  analytical  estimate  of  the  viscous  effects  of  hypersonic  flow 
over  an  insulated  wedge.  Shen  assumed  that  the  viscous  effects 
are  confined  in  tne  region  from  the  body  surface  to  the  leading 
edge  shock  wave.  He  concluded  [2]  that  the  compressible 
boundary  layer  equations  are  valid,  as  a  first  approximation,  in 
the  viscous  layer  on  the  body.  By  a  simple  momentum  integral 
analysis,  it  was  possible  to  treat  the  shock  wave  and  boundary 
layer  interaction  phenomena  at  hypersonic  speeds  as  a  complete 
problem.  In  an  experimental  study  of  the  hypersonic  flow  on  a 
flat  plate  at  the  Langley  NACA  Laboratory,  probably  concurrent 
to  Shen’s  research,  Becker  [3]  observed  that,  at  M*  =  6.86, 
there  is  a  significant  pressure  rise  on  the  frontal  portion  of  the 
plate.  Bertram  [4]  in  an  attempt  to  interpret  these  experimen¬ 
tal  findings,  proposed  in  1952  an  approximate  method  for  de¬ 
termining  the  displacement  effects  and  viscous  drag  of  laminar 
boundary  layers  in  two-dimensional  hypersonic  flow.  In 
Bertram’s  work,  the  region  behind  the  shock  wave  was  treated 
as  consisting  of  a  viscous  boundary  layer  and  an  inviscid  region 
between  the  shock  and  the  boundary  layer.  In  1952,  Lees  and 
Probstein  [5]  at  Princeton  gave  a  more  rigorous  treatment  of 
Bertram’s  problem  of  weak  interaction  between  the  shock  wave 
and  the  boundary  layer  in  hypersonic  flow.  In  the  same  year, 
the  present  authors  presented  a  revised  theory  [6]  of  Shen’s 
problem  of  strong  interaction.  In  [6],  it  was  discovered  that 
(1)  the  strong  interaction  problem  yields  a  pressure  distribution 

law  in  the  front  part  of  a  flat  plate,  pctx"^  ,  (2)  the  boundary 
layer  in  the  strong  interaction  region  grows  like  dax3/4,  and 
(3)  the  solutions  have  the  characteristics  of  “similar”  solutions 
of  the  compressible  boundary  layer  equations.  These  findings 
were  later  confirmed  by  Lees  [7],  who  formulated  an  analytical 
theory,  not  using  the  momentum  integral  method,  of  the  strong 
interaction  problem.  However,  in  [7],  Lees  retained  the 
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assumption  that  the  shock  wave  and  the  boundary  layer  are 
separated  by  an  inviscid  rotational  flow  region.  In  order  to  af¬ 
fect  a  reasonable  matching  between  the  viscous  and  inviscid 
regions,  Lees  had  to  introduce  the  tangent-wedge  approxima¬ 
tion  [8] . 

In  an  entirely  general  discussion  of  a  method  of  treating  the 
compressible  boundary  layer  equation,  Stewartson  [9]  pointed 
out  in  1949  the  existence  of  “similar”  solutions,  Stewartson’s 
concept  was  carried  out  in  great  detail  by  the  present  authors 
in  1952  [10].  A  new  approach  to  the  problem  of  “similar”  solu¬ 
tions  was  also  presented  in  a  recent  paper  [11]  in  which  it  was 
pointed  out  that  Lees’  results  in  (7)  are  a  special  case  of  the 
“similar”  solutions.  Thus,  the  treatment  of  the  hypersonic  in¬ 
teraction  problem  and  the  discussion  on  “similar”  solutions 
are  completely  correlated. 

It  is  interesting  to  note  that  the  hypersonic  interaction  phe¬ 
nomena  give  a  physically  significant  example  of  the  “similar” 
solutions  of  the  compressible  boundary  layer  equations,  but  the 
most  interesting  point  is  that  by  the  use  of  the  “similar”  solu¬ 
tions  it  is  possible  to  extend  the  analytical  treatment  of  the  hy¬ 
personic  interaction  problem  to  include  more  general  cases. 
For  instance,  for  the  strong  interaction  problem,  there  is  not 
yet  a  satisfactory  theory  dealing  with  the  case  of  a  noninsulated 
flat  plate.  This  particular  problem  will  be  treated  in  the  pres¬ 
ent  paper  by  the  method  of  “similar”  solutions. 

The  adoption  of  the  present  treatment,  however,  depends  on 
two  vital  conditions:  (1)  In  the  case  M  «,  >  >  1,  the  velocity  at 
the  edge  of  the  boundary  layer  is  nearly  constant,  ux  —  u^; 

(2)  In  the  frontal  part  of  the  flat  plate,  pc*x-1/2.  Recent  experi¬ 
mental  studies  at  GALCIT  [12]  and  Princeton  [13]  have  con¬ 
firmed  these  conditions.  They  can  be  readily  utilized  to  reduce 
the  system  of  partial  differential  equations  to  ordinary  equa¬ 
tions  by  the  method  of  [11],  Thus,  a  theory  of  hypersonic  vis¬ 
cous  flow  on  a  noninsulated  flat  plate  can  be  presented  entirely 
within  the  framework  of  the  method  of  “similar”  solutions. 

The  results  of  the  present  study  are  of  practical  importance 
because  they  give  an  estimate  of  the  effects  of  heat  transfer  in 
the  strong  interaction  region  on  a  flat  plate.  Presumably,  if 
sufficient  knowledge  of  the  heat  transfer  effects  is  accumulated, 
it  will  be  possible  eventually  to  control  effectively  the  strong 
interaction  effects. 
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In  f  1 1  ],  it  has  been  shown  that  if  (1)  the  Prandtl  number  of 
the  fluid  is  assumed  unity,  and  (2)  the  viscosity  of  the  fluid  is 

U  Q 

assumed  to  behave  like  a  linear  function  of  T,  i.e.,  fi  =  — -  T, 

1  o 

then  the  partial  differential  equations  of  the  steady,  two- 
dimensional  compressible  boundary  layer  can  be  transformed 
into  a  system  of  ordinary  equations  as  follows: 

KK"  +  K'"  =  —  (K'2  -  G)  (1) 

m 

G"  +  KG'  =0  (2) 


provided  that  the  main  flow  outside  the  boundary  layer  satisfies 
the  following  differential  equation: 


In  Eqs  1  and  2,  K  and  G  are  functions  of  a  similar  variable: 


ri  (x,  y)  = 


1-m  dy 

1  T(x7y) 


(4) 


and  the  prime  denotes  differentiation  with  respect  to  r) .  The 
physical  significance  of  K  and  G  is  evident  from  the  following 
definitions: 


K'  =  — y— r  ,  G  =  £ 

ux(x)  Hj 

The  proper  boundary  conditions  are 

ri  =  o,  K  =  K1  =  o,  G  =  —■ 

1  o 

T)  >oo,  K 1  =  1 ,  G  =  1 


(5) 


(6) 

(7) 


In  hypersonic  flow,  >  >  1,  the  flow  velocity  behind  the 
leading  edge  shock  wave  is  not  much  different  from  the  undis¬ 
turbed  free  stream  velocity,  i.e., 
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For  the  “strong*  interaction  problem  of  the  hypersonic  viscous 
flow  on  a  flat  plate,  if  the  outer  edge  of  the  boundary  layer  is 
regarded,  to  a  first  approximation,  as  a  streamline,  then  it  can 
be  shown  [14]  that  Eq  3  becomes 


It  is  concluded,  therefore,  that  for  the  hypersonic  boundary 
layer  problem,  the  constant  parameter  m  in  Eq  1  must  satisfy 
the  following  condition: 

m -- (t^y)  i'-r1)  (1° 

where  n  is  an  index  number  of  the  main  stream  pressure,  viz., 

^  a  xn.  Previous  calculations  in  [6]  show  that  n  =  -^ 

Po  “ 

Hence 


Eqs  1  and  2  with  the  boundary  conditions  of  Eqs  6  and  7,  in  the 
case  of  values  for  m  as  given  in  Eq  11,  have  been  integrated  on 
a  REAC  machine  with  Y  =  1.4  and  Y  =  1.67  for  a  range  of 
values  of  Tw/T0.  These  solutions  have  previously  been  re¬ 
ported  in  [11],  In  the  present  paper,  this  data  will  be  used  to 
compute,  for  various  Tw/T0’s,  the  hypersonic  viscous  flow 
characteristics  on  a  flat  plate  in  an  air  atmosphere  or  helium 
atmosphere.  These  results  are  invaluable  in  the  estimation  of 
the  effectiveness  of  surface  cooling  or  heating  in  controlling  the 
“strong"  interaction  effects. 

Computations  in  [14]  show  that  within  the  framework  of  an 
approximation  by  the  method  of  similar  solutions,  the  “strong" 
interaction  between  the  leading  edge  shock  wave  and  the  viscous 
boundary  layer  on  a  flat  plate  in  hypersonic  flow  is  character¬ 
ized  by  the  following  general  laws: 
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VHe. 


'  ml  Y/a 

y*rj 
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(12) 

(13) 


CfcoV^  = 


(14) 


ch  „oV^ex 


(15) 


where  A,  B,  C,  D  are  defined  as  follows: 


The  values  of  A,  B,  C,  D  for  various  Tw/T0’s,  with  Y  =  1.4  and 
Y  =  1.67,  are  presented  in  Figs.  1,  2,  3,  and  4,  respectively. 
They  are  also  tabulated  in  Table  I.  It  is  observed  that  the 
integral 


I  = 


L 


(G  -  K'2)d 77 


(20) 


280 


CONFERENCE  ON  FLUID  MECHANICS 
Table  I 


VALUES  OF  A,  B,  C,  D  FOR  VARIOUS  Tw/T0’s 


Tw/To 

y  =  1.4 

A  B  C  D 

A 

y  =  1.67 

BCD 

0 

.1485  .3967  .2082  .1888 

.2614 

.4572 

.2867  .2494 

.2 

.2316  .4954  .2793 

.4029 

.5675 

.3946 

.6 

.3773  .6323  .4118  .3105 

.6596 

.7261 

.6015  .4199 

1.0 

.5140  .7381  .5488  .3732 

.9207 

.8577 

.8199  .5042 

2.0 

1.506 

1.097 

1.361  .6832 

enters  all  the  formulae  for  A,  B,  C,  and  D.  This  integral  is 
evaluated  by  Simpson’s  rule,  with  At]  =  0.1,  from  the  data  ob¬ 
tained  from  the  REAC.  The  values  of  I  for  various  Tw/T0’s 
are  given  in  Table  II. 


VALUES  OF  I 


Table  II 

c© 

(G  -  K,2)dr?  FOR  VARIOUS  Tw/T0’s 


y 

Tw/T0 

I 

1.4 

0 

.3820 

1.4 

.2 

.5956 

1.4 

.6 

.9703 

1.4 

1.0 

1.322 

1.67 

0 

.3507 

1.67 

.2 

.5404 

1.67 

.6 

.8848 

1.67 

1.0 

1.235 

1.67 

In  Figs.  1-4,  Tw/Tq  =  1.0  corresponds  to  an  insulated  sur¬ 
face;  Tw/T0  <  1.0,  a  cooled  surface;  and  Tw/T0  >  1.0,  a 
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.  Values  of  C  (Eq.  18)  vs.  Tw/T0,  y  «  1.4  and  y  =  1.67 
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.  Values  of  D  (Eq.  19)  vs.  Tw/T0,  y  =  1.4  and  y  =  1.67 
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heated  surface.  From  these  figures,  it  can  readily  be  concluded 
that  cooling  the  surface  will  decrease  the  strong  interaction  ef¬ 
fects  between  the  shock  wave  and  the  boundary  layer.  The  ef¬ 
fects  of  cooling  are  many  fold:  (1)  cooling  the  surface  tends  to 
thin  down  the  boundary  layer  thickness,  (Fig.  2);  (2)  cooling  the 
surface  tends  to  induce  a  smaller  pressure  rise  in  the  inter- 
a-.t  .on  zone  (Fig.  1);  (3)  cooling  the  surface  tends  to  lower  the 
>kin  friction  coefficient  for  a  fixed  Rex  (Fig.  3).  The  values  of 
2Choo/Cfoo  for  various  Tw/T0’s  are  presented  in  Fig.  5.  This 

parameter  tends  to  increase  when  Tw/T0  decreases  to  zero. 

The  effects  of  heating  the  surface  are  just  the  opposite  of  the 
above. 


Discussion 

The  recent  experimental  findings  at  GALCIT  [12]  confirm 
that 

(1)  the  static  pressure  on  the  plate  surface  varies  as  x"1/* 
(Fig.  6) 
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6.  Static  Pressure  Distribution  on  a  Flat  Plate. 


(2)  the  velocity  at  the  edge  of  the  boundary  layer  is  nearly 
constant:  u,  :  ux 

On  the  basis  of  these  facts,  a  “similar*  solution  can  be  found 
for  the  problem  erf  the  hypersonic  strong  interaction  between  the 
shock  wave  and  the  boundary  layer  on  a  flat  plate.  More  ex¬ 
tensive  experimental  studies  are  necessary  to  verify  the  theo¬ 
retical  predictions  of  the  behavior  of  the  boundary  layer  skin 
friction  and  heat  transfer  characteristics. 

At  the  present  time,  in  the  particular  case  of  an  insulated 
flat  plate,  a  comparison  of  the  values  of  A  from  experiment  and 
from  theory  can  be  made.  This  comparison  is  shown  in  Fig.  7. 
It  is  seen  that  the  air  data  at  =  5.8  obtained  in  the  GALC1T 
Hypersonic  Wind  Tunnel,  Leg  No.  1,  seem  to  agree  quite  well 
with  the  theories.  It  is  also  interesting  to  note  that  the  agree¬ 
ment  with  the  theory  in  [6]  seems  even  better  than  the  present 
computations.  This  is  perhaps  due  to  two  inaccuracies  in  the 
present  analysis: 

(1)  the  assumed  linear  variation  of  \i  with  T, 

(2)  the  approximate  matching  of  the  viscous  and  inviscid  re¬ 
gions  by  the  tangent-wedge  approximation. 
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7.  Comparison  of  Theories  and  Experiments. 


In  [6],  of  course,  the  more  accurate  variation  of  fi  with  T  is 
used,  namely,  fi  a.  Tw,  with  the  result  that  the  hypersonic 

2+ Co  I  |  * 

parameter,  '\Rex »  *s  obtained  instead  of  the  more  re¬ 
strictive  parameter,  M*^R^  .  Stewartson,  in  a  recent  study 

[15]  of  the  hypersonic  strong  interaction  phenomena,  assumed 
that  a  “similar"  solution  also  exists  for  the  inviscid  rotational 
flow  region  between  the  shock  wave  and  the  boundary  layer. 
Thus,  he  was  able  to  affect  a  matching  at  the  edge  of  the  bound¬ 
ary  layer  without  using  the  tangent-wedge  approximation.  With 
Stewartson’s  method  of  matching,  the  values  of  A  would  come 
closer  to  the  experimental  values. 

The  helium  data  at  Moo  =  12.7  obtained  in  the  Princeton 
tunnel  [13]  give  considerably  larger  values  of  A  than  the 
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theoretical  predictions,  which  might  be  due  to  the  fact  that  the 
flat  plate  model  used  at  Princeton  has  a  larger  leading  edge 
radius.  The  importance  of  the  effects  of  the  bluntness  of  the 
leading  edge  has  been  recently  brought  out  by  Bertram  [16]  and 
Lees  [17].  As  pointed  out  in  [17],  the  present  "similar*  solu¬ 
tions  probably  represent  the  correct  limiting  forms  of  the  solu¬ 
tion  for  the  case  of  an  ideally  sharp  leading  edge  plate  when  the 
hypersonic  parameter,  M%0/^Rex,  is  held  constant  as  +  <». 

In  this  idealized  situation,  the  consideration  of  the  bluntness  of 
the  leading  edge  and  the  estimation  of  the  effects  of  the  entropy 
layer  from  the  leading  edge  region  are  not  possible.  To  remedy 
this  unrealistic  omission,  it  is  necessary  to  study  the  leading 
edge  region  more  carefully. 

At  the  present,  nothing  can  be  said  about  the  comparison  of 
experiment  and  theory  for  the  case  of  hypersonic  strong  inter¬ 
action  on  a  noninsulated  flat  plate.  Extensive  experimental 
studies  in  this  direction  should  be  initiated.  Nevertheless,  by 
inference  from  the  present  results,  it  is  felt  that  the  hyper¬ 
sonic  strong  interaction  phenomena  can  be  reasonably  con¬ 
trolled  by  providing  sufficient  surface  cooling. 
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